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Введение. Цель работы — повышение надёжности методов 
прогнозирования поведения гидроксида лития в пароводя- 
ном тракте на ТЭС и АЭС и оперативного контроля водо- 
родного показателя паровых растворов. Разработан метод 
оперативного контроля на основе кондуктометрических 
измерений водородного показателя паровых растворов ин- 
гибитора коррозии конструкционных материалов гидрокси- 
да лития, применяемого на ТЭС и АЭС. 

Материалы и методы. Для практической реализации высо- 
котемпературного оперативного контроля водородного по- 
казателя паровых растворов использован метод математиче- 
ского моделирования. 

Результаты исследования. Разработан метод контроля во- 
дородного показателя паровых растворов гидроксида лития, 
основанный на определении концентрации в паре путём 
конденсации пара в охлаждаемом кондуктометрическом 
датчике, размещённом в паровом пространстве парогенера- 
тора, что позволило значительно повысить точность опреде- 
ления концентрации гидроксида лития. Предложены урав- 
нения, описывающие изменение предельной молярной экви- 
валентной электропроводности и констант диссоциации 
гидроксида лития в широком диапазоне изменения парамет- 
ров состояния. 

Обсуждение и заключения. Предложенные метод оператив- 
ного контроля водородного показателя паровых растворов и 
математические модели для расчёта предельной эквивалент- 
ной молярной электропроводности и констант диссоциации 
обеспечивают приемлемую для практических расчётов по- 
грешность и возможность автоматизации измерений. При 
увеличении температуры пара до 573,15 К возникает необ- 
ходимость в увеличении концентрации гидроксида лития в 
паре до 10? молыкг. 


Ключевые слова: машиностроение, коррозионные процес- 
сы, автоматический контроль, паровая фаза, кондуктометри- 
ческие методы, мониторинг. 
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Введение. При организации и ведении водно-химических режимов (ВХР) на ТЭС и АЭС одной из 
главных задач является сведение к минимуму коррозионных повреждений теплосилового оборудования в зоне 
контакта с водой и паром [1]. Дозирование в питательную воду [4ОН в необходимых количествах обеспечивает 
создание на поверхности металла прочной защитной литий-ферритной плёнки, замедляющей коррозию и кор- 
розионное растрескивание теплообменных труб парогенераторов АЭС с ВВЭР-1000 [2, 3], а также значитель- 
ное увеличение времени наработки оборудования на отказ [3]. При дозировании МОН в питательную воду ба- 
рабанных котлов ТГ-104 на ТЭС также наблюдается формирование защитной литий-ферритной плёнки [4]. На 
АЭС ЦОН применяют в России и за рубежом [4—6]. В перспективных атомных реакторах ГУ поколения $5СИК 
(бирегстиса ижиет-сое4 геас1о’), обеспечивающих увеличение термического КПД по сравнению с существу- 
ющими (РР) с 33 до 44 %, при давлении Р=25 МПа температуру Т водяного пара планируют увеличить до 
953 К [6]. Для создания на ТЭС и АЭС оперативного контроля значений рН водных и паровых растворов [4ОН, 
определяющих эффективность защиты от коррозии, необходимы данные о предельной эквивалентной электро- 
проводности Л. и константах диссоциации К/ 4ОН при высоких Р и Т, представленные в виде, удобном для 
практических расчётов. Методика расчёта рН растворов Г[4ОН на основании данных об удельной электропро- 
водности Хх, предложенная авторами для жидкой фазы [7], неприемлема для паровой фазы при малых значениях 
плотности р и концентрации С в связи с высокой погрешностью определения хи К/ [8]. 

Целью настоящей работы является повышение надёжности методов прогнозирования поведения ОН 
в пароводяном тракте на ТЭС и АЭС и оперативного контроля рН паровых растворов. 

Для достижения поставленной цели созданы уравнения для расчёта Ло и Ки на линии насыщения и в 
перегретом паре при Т до 773,15 К, а также разработана методика оперативного контроля рН паровых раство- 
ров [4ОН. 

Материалы и методы. По сравнению с традиционным контролем, связанным с отбором и охлаждени- 
ем проб, оперативный контроль рН при помощи кондуктометрических датчиков, установленных в теплоноси- 
теле парогенераторов и работающих при Ти Р, близких к рабочим, является более объективным и практически 
безынерционным [7, 8]. В расчётах рН для жидкой фазы в [7] были использованы экспериментальные данные о 
Ка ЧОН [8], ионном произведении К, [9], и диэлектрической проницаемости & для воды [10]. Определение кон- 
центрации С водных растворов МОН осуществлялось на основании зависимостей С МОНЕ} (х МОН)т-сотзь ПО- 
лученных на основании экспериментальных данных о // растворов. 

При измерении Хх растворов в паровой фазе (особенно при малых ри С), когда полярные свойства рас- 
творителя ослаблены, степень диссоциации 4ОН мала и сопротивление межэлектродного пространства изме- 
рительной ячейки кондуктометрического преобразователя может стать соизмеримым с сопротивлением элек- 
трического изолятора, что приведёт к большим погрешностям определения Хх [8]. 

Для практической реализации высокотемпературного оперативного контроля рН необходимые данные 
о Льи К; могут быть получены в исследуемом диапазоне изменения параметров состояния методом математи- 
ческого моделирования и представлены в виде уравнений. 

Результаты исследования. На основании анализа существующих методов в данной работе предложен 
метод контроля рН паровых растворов /[4ОН, основанный на определении С ШОН в паре путём конденсации 
пара в кондуктометрическом датчике с охлаждаемым капилляром [11], размещённом в паровом пространстве 
парогенератора в потоке теплоносителя. В этом случае С МОН в паре равна С Г4ОН в конденсате пара и может 
быть определена из зависимостей С МОН=(х МОН)т-сопя, полученных для жидкой фазы [7]. Конденсация пара 
позволяет значительно повысить точность определения х и С [8]. Вычисления рН для парорастворов [4ОН про- 
изводились по методике, описанной в [12], с использованием данных [8] о Ки. 

С целью оценки эффективности применения /4ОН для защиты от коррозии в атомных реакторах ТУ по- 
коления 5СУ авторами [6] на основании экспериментальных данных о К; [13-15] выполнен расчёт рН водных 
и паровых растворов [4ОН при Р=25 МПа в широком диапазоне изменения р, Т, С и предложено эмпирическое 
уравнение для вычисления моляльной константы ассоциации К, [16] в диапазоне изменения р от 50 до 1000 
кг/м” и Гот 373 до 873 К: 

08 1оКа = а, + аз + + (аз + а4/т)-108 106 + (@5 + ав). (1085), 
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где б=р/рс; т=Т/То; ре, Те — плотность и температура раствора в критической точке; а1, а>, аз, ал, а5, ав — кон- 
станты. 

В перегретом паре с уменьшением ри ТГ К, ЧОН уменьшается, уменьшается и разность между рН Н›О 
и рН парорастворов Г[4ОН, что снижает эффективность защиты от коррозии элементов конструкций при высо- 
ких Ри Г. Авторами вычислены значения рН парорастворов [4ОН при Р=25Мпа; Т=673, 723 и 773 К для С=10°, 
105, 10“, 103 и 10 ? молыкг. В расчётах были использованы экспериментальные данные о Ка [4ОН [8], ионном 
произведении К,, [9], и диэлектрической проницаемости воды [10]. На рис. 1 представлена зависимость рН 
водных и паровых растворов [4ОН от р при Р=25 МПа для С=10`°-10?, молыКг. 

В области р<330 кг/м” авторы [13—15] не выполняли экспериментальные исследования, поэтому на ри- 
сунке представлены только результаты расчёта рН ГОН с использованием данных о К, [4ОН, полученных по 
уравнению авторов [16 ]. 

773 и 9 
723 
12.0 873] 673 623 573 523473 373 











5,0 


100 200 300 400 500 600 700 800 900 1000 
р, кг/мз 


Рис.1. Зависимость рН водных и паровых растворов /ЛОН от плотности р при Р=25 МПа: 
1— Н.0; 2, 3, 4, 5,6 — С МОН=10°, 10°, 10%, 10°, 10? молыкг; 
А,о — расчёт рН по данным о Ки [13—15] и [8]; 
линии 2, 3, 4, 5, 6 — расчёт рН по данным о К, полученным из уравнения [16] 


Полученные авторами ранее на основании измерений электропроводности высокотемпературных рас- 
творов /[4ОН, значения Ло и К. [8], необходимые для создания уравнений, описывающих изменение Л. и К. при 
изменении р и Т, приведены в таблицах | и 2 в сопоставлении с данными других авторов. Для создания уравне- 
ний, описывающих поведение Л. на линии насыщения, с целью уменьшения погрешности аппроксимации дан- 
ных диапазон изменения р был разбит на 2 части (от 100 до 500 и от 500 до 958 кг/м), для каждой из которых 
предложено уравнение вида: 

Ло=Бо:р* + В1-р* + Бур + Бь, 
где Бе, Би, Б›, Бз — определяемые константы; р — плотность раствора, кг/м?. Значения Бо, Ба, Б>, Бз для р=(100- 
500) и (500-958) кг/м* равны -—0,000004; 0,005878; —2,810820; 1942,943226 и —1,2571 1-10; 0,024177411; — 
16,58559801; 5334,222549 соответственно. Для значений рК.=ю0®1 Ка вычислены значения коэффициентов 
4-4 в уравнении, аналогичном представленному в [16] и приведённому выше: 
рК=а, + ат + (4 + а/®)1овтоб + (45 + а/®)- (1035), 

где б=р/рс; т=Т/Т,; ре, Т. — плотность и температура раствора в критической точке; 41, 4, 43, 44, 45, 46 — кон- 
станты. 

При вычислении констант 4-4 были использованы полученные авторами ранее экспериментальные 
данные о К/ [8]. На линии насыщения авторы учитывали только свои экспериментальные данные. Расчёты для 
перегретого пара при Т=673,15 и 773,15 К выполнены в двух вариантах: в первом рассматривались только ав- 


торские данные о К/ [8], во втором — все данные, приведённые в табл. 2. 
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Таблица 1 


Зависимость предельной молярной электропроводности Л, 4ОН 
от плотности р на линии насыщения и в перегретом паре 

































































р, кг/м? Ло.10, См-м”/кмоль 
Линия насыщения Т= 673 К Т=773К 
100 17109[8] 1675[8] 1665[8] 
130 1675[8] 1630[8] 1620[8] 
160 1630[8] 1590[8] 1575[8] 
200 1575[8] 1550[8] 1545[8] 
300 1510[8] 1480[8] 1475[8] 
400 1475[8] 1460[8]; 1392[13] 1455[8] 
433 1452[15] 
500 1450[8] 1435[8]; 1539[13] 1428[8]; 1574[13] 
600 1370[8] 1380[8]; 135513] 1388[8]; 137713] 
621 127715] (Ф = 23,86 МПа) - - 
700 1265[8] 1275[8]; 127013] 1276[13] 
714 1189[15] (Р=9,52 МПа) - - 
800 1100[8] 1088[13] 1109[13] 
871 1017[15] (Р = 9,91МПа) - - 
917 78018] - - 
922 754[15] (Р =9,4 МПа) - - 
950 594[13] - - 
958 573[8] - 

















Для линии насыщения: 4,=119,50984; 4=-115,333 


4=183,09356. 
Для первого варианта расчёта для перегретого пара: 4\=7,6323975; 4=3,214114; 4-=10,9289; 


4=2,5463453; 45=-10,34132; ф=15,190972; для второго варианта расчёта 4=10,113404; 4›=—6,133638; 4= 


17,85984; 4.=8,888042: 4.=-8,405401:; 4=14,648395. 
Отклонения в % экспериментальных значений Л. и рК. от вычисленных по предложенным уравнениям 
представлены в табл. 3. 


47; 4з=174,20233; 4=167,90533; 45=—143,78368; 










































































Таблица 2 
Зависимость отрицательного логарифма констант диссоциации [4ОН рКа 
от плотности р на линии насыщения и в перегретом паре 
3 РКа, моль кг 
р, кг/м 
Линия насыщения Т = 673,15 К Т= 773,15 К 

100 10,04[8] 10,07[8]; 11,06 [16] 10,16[8]; 13,34 [16] 
130 8,4118] 8,56[8] 8,67[8] 
160 7,24 [8] 7,32[8] 7,69[8] 
200 6,25[8] 6,42[8]; 6,45[16] 6,89[8]; 8,16[16] 
300 4,40[8] 4,69[8]; 4,13[16] 5,20[8]; 5,64[16] 
379 3,4815] 
400 3,86[8] 3,98[8]; 2,82[16]; 3,14[13] 4,3118]; 4,0116] 
433 3,1115] 
500 3,28[8] 3,38[8]}; 2,8113] 3,4318]; 3,15[13] 
600 2,5218]; 2,45[8]; 2,32[13]; 1,26[16] 2,58[8]; 2,43[13]; 1,97[16] 
650 2,13[13] 2,25[13] 
700 2,08[8]; 2,0118]; 1,83[13] 2,07[13] 
712 1,59[17] 
750 1,82[13] 1,89[13] 
800 1,63[131; 0,343 [16] 1,65[131; 0,7536] 
807 1,4115] (Р= 1.02 МПа) - - 
922 1,29[15] (Р=9,4 МПа) - - 
961 1,0315] (Р=4,88 МПа) - - 
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Отклонение (в %) экспериментальных значений предельной эквивалентной 





Таблица 3 


электропроводности (6Л.,) и отрицательного логарифма моляльных констант диссоциации (брКа) 


от значений, вычисленных по уравнениям 
























































Л. брКа 
р Т=673,15 К Т=773,15 К Т=673,15 К Т=773,15 К 
ди Линия Линия 
кг/м 1 вариант ]вариант 2 вариант 2 вариант 
насыщения насыщения в ы а ы 
расчёта расчёта расчёта расчёта 
3,40[8]; 10,981; 
100 0,36[8] 0,0418] 1,208] 1,36[8] 
5,86 [16] 15,5[16] 
130 0,47[8] 0,23[8] 0,18[8] 1,22[8] 2,42[8] 11,78] 
160 0,28[8] 0,3318] 2,62[8] 1,86[8] 3,42 [8] 10,5[8] 
0,52[8]1; 5,95[8]; 
200 0,34[8] 0,1118] 0,92[8] 0,39[8] 
0,98[16] 10,5[16] 
3,86[8]; 1,9981; 
300 0,09[8] 2,50[8] 0,29[8] 
9,17[16] 5,96[16] 
379 3,05[15] 
14,8[8]; 
6,44[8]; 
400 0,0318] 0,12[8] 4,99[8] 4,07[8] 20,2[16]; 
0,56[16]; 
7,96[13] 
433 0,4115] 
21,778]; 8,78[8]; 
500 0,09[8] 0,16[8] 8,95[8] 2,40[8] 
5,90[13] 0,67[13] 
13,7218]; 5,24[8]; 
600 0,10[8] 0,17[8] 4,56[8] 6,1118] 8,89[131; 0,60[13]; 
67,8[16 24,1[16] 
650 - - - - 10,6[13] 4,02[13] 
14,38[8]; 
700 0,44[8] 0,05[8] 7,828] - 8,0313] 
5,95[13] 
750 - - - 14,213] 11,5[13] 
800 0,26[8] - - - 12,8[13] 11,3[13] 
917 2,27[8] - - - - - 
958 1,5118] - - - - - 





























Обсуждение и заключения. Как очевидно из рис. 1, полученные авторами значения рН 4ОН при 
Т=673 К ниже по сравнению с полученными путём расчёта [6] в среднем на 0,2; при Т=723 К выше на 0,1-0,2; 
при 7=773 К выше на 0,1-0,7 единиц рН. 

Таким образом, величина рассогласования данных о рН ШОН в перегретом паре, полученных авторами 
и изложенных в [6], возрастает с увеличением Т и С. Такое рассогласование можно считать удовлетворитель- 
ным, если учесть, что с уменьшением р до 200-100 кг/м* погрешность экспериментальных данных [8] о Ка 4ОН 
возрастает до 180 % и, кроме того, при составлении уравнения авторами [6] не использовались эксперимен- 
тальные данные о Ки при р<330 кг/м‘. 

Отклонение экспериментальных значений Л. и рК/ от вычисленных по предложенным в настоящей ра- 
боте уравнениям на линии насыщения не превышает 1 % и только при р=917 и р=958 кг/м* достигает 2,27 и 
1,51 % для Л. (табл. 3). В перегретом паре при 7=673,15 и Т=773,15 К для первого варианта расчёта с использо- 
ванием только авторских экспериментальных данных о К/ [8] наибольшее отклонение экспериментальных дан- 
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ных рКо от расчётных не превышает 2,5 % при р=100-300 кг/м? и возрастает до 8,95 % с увеличением р до 400- 
700 кг/м. Для второго варианта расчёта с учётом данных всех авторов, представленных в табл. 2, наибольшая 
погрешность расчёта по уравнению при р=100-300 кг/м и 7=673,15 К; Т=773,15 К достигает 15,5 % для данных 
авторов [16], а при р=100-300 кг/м* и Т=673,15 К; Т=773,15 К — 67,8 %. 

На основании анализа полученных результатов можно сделать вывод о том, что использование ОН в 
качестве корректирующей добавки для регулирования рН с целью предотвращения коррозионного разрушения 
элементов конструкций перспективных атомных реакторов ГУ поколения 5СУЛ не вызывает затруднений при 
7Т<673,15 К и Р=25 МПа. С ростом температуры уже при Т=773,15 К для эффективной защиты от коррозии 
необходимо увеличивать С [4ОН до 107 моль/кг, что выдвигает задачу исследования растворимости 4ОН в 
паре при этих параметрах состояния. 

Предложенные в настоящей работе метод оперативного контроля рН паровых растворов ОН и мате- 
матические модели для расчёта Л. и Клобеспечивают возможность определения рН при помощи измерительно- 
вычислительного комплекса, работающего непрерывно в автоматическом режиме, имеющего малое время за- 
паздывания сигнала и обеспечивающего приемлемую для практических расчётов погрешность результатов из- 
мерения. Использование таких устройств в современных системах химико-технологического мониторинга на 
ТЭС и АЭС позволит значительно повысить надёжность методов прогнозирования поведения /4ОН в парово- 
дяном тракте и оперативного контроля рН. 
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